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We discuss a new way to develop the exactly solvable model of the QCD critical endpoint by 
matching the deconfinement phase transition line for the system of quark-gluon bags with the line 
of their vanishing surface tension coefficient. In contrast to all previous findings in such models the 
deconfined phase is defined not by an essential singularity of the isobaric partition function, but by 
its simple pole. As a result we find out that the first order deconfinement phase transition which is 
defined by a discontinuity of the first derivative of system pressure is generated by a discontinuity 
of the derivative of surface tension coefficient. 
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I. INTRODUCTION 

Recently intensive theoretical and experimental search 
for the (tri-) critical endpoint of strongly interacting mat- 
ter at small enough chemical potential become very fas- 
cinating and promising branches of research activity in 
the context of relativistic heavy ion programs of many 
laboratories. In particular, the most powerful computers 
and very sophisticated algorithms are used for the lattice 
quantum chromodynamics (LQCD) simulations to locate 
this endpoint with maximal accuracy and to study its ori- 
gin and properties [1, 2, 3 , but despite these efforts the 
present results are still far from being conclusive. The 
general arguments for a similarity with the critical point 
features in the other substances are not much convincing 
mainly because of a lack of rigorous critical point theory 
which exists, in a sense, for the spin systems only [4^, 
whereas the origin and physics of the critical point for 
realistic gases and nuclear matter are described, at best, 
phenomenologically. For example, the Fisher droplet 
model (FDM) O [6] turns out rather efficient in studying 
the critical point of realistic gases. This model was ap- 
plied to many different systems with the different extents 
of success including a nuclear multifragmentation [7 , a 
nucleation of real fluids [8 , the compressibility factor of 
real fluids [9 , the clusters of the Ising model [IQHH] and 
the percolation clusters [12 , but really its phase diagram 
does not include the fluid at all and, therefore, is not 
completely satisfactory and theoretically well defined. 

The statistical multifragmentation model (SMM) [13] 
looks much more elaborated in this aspect because de- 
fines the phase diagram of the nuclear liquid- vapor type 
phase transition (PT) in some controlled approximation 
[m [15] and predicts the critical (tri-critical) endpoint ex- 
istence for the Fisher exponent 0<r<l(l<r<2) [T5] 
together with giving the possibility to calculate the corre- 
sponding critical exponents [16 . However, the predicted 
location of the SMM (tri-)critical endpoint at maximal 
density of the nuclear liquid does not seem to be quite 
realistic. Actually, the relations between all these criti- 
cal points are not well established yet. In principle, the 
Complement method [TP provides us with the possibil- 



ity to describe accurately the size distribution of large 
clusters of 2- and 3-dimensional Ising model within the 
FDM in rather wide temperature interval but the de- 
tailed numerical comparison of the Ising model and the 
FDM critical endpoints is hardly possible because of the 
large fluctuations even in relatively small systems. In 
meantime, taking the formal limit of the vanishing nu- 
cleon proper volume leads to the situation in which the 
SMM grand canonical partition function covers the FDM 
partition function, but the analytical properties are not 
the same and, as a result, the condensation particle den- 
sity of gaseous phase is finite not for the Fisher exponent 
r < 2, as in the SMM, but for r > 2 [16 . Besides, it 
leads to the various correlations between the r exponent 
and other critical indices in the FDM [5 and the corre- 
sponding relations in the SMM are different what signals 
the universality classes for these models are different as 
well [16 . 

As to the model calculations of the QCD phase struc- 
ture they are based on the universality arguments ad- 
vanced in [18] and concern mainly the temperature driven 
chiral symmetry restoration transition. In fact, it can 
not provide the reliable conclusion about the transition 
order at /i = 0, its dependence on the number of flavors 
[19 and especially about the location of the point (tri- 
critical) on the PT line where the transition changes its 
order. It seems the lattice QCD (LQCD) simulations at 
vanishing /i give more definite evidence that this temper- 
ature driven phenomenon could really be a crossover [17] . 
Then, clearly the (/i, T) phase diagram contains a critical 
point caused by the /i = crossover turning into the first 
order PT. However, this wisdom is rather questionable 
as well [20 . 

Furthermore, the recent LQCD simulations [![ [21 [2T] 
teach us that even at high temperatures up to a few Tc 
(Tc is the cross-over temperature), a QGP does not con- 
sist of the weakly interacting quarks and gluons and its 
pressure and energy density are well below of the cor- 
responding quantities of non-interacting quarks and glu- 
ons. Although such a strongly coupled QGP (sQGP) [22] 
has put a new framework for the QCD phenomenology, 
the feasibility of understanding such a behavior within 



2 



the AdS/CFT duahty[23^ or statistical approaches is far 
from being simple and transparent. 

Here we investigate the possibility to resolve the prob- 
lem by formulating an approach based on the model 
of quark-gluon bags with surface tension (QGBSTM) 
[24l [25l [26] . The paper is organized as follows. Sect. II 
contains the formulation of model basic elements (here- 
after the model is named QGBSTM2 in order to distin- 
guish it from the model with the tri-critical endpoint). In 
Sect. Ill we analyze all possible singularities of the QGB- 
STM2 isobaric partition for non- vanishing baryonic den- 
sities and discuss the necessary conditions for the criti- 
cal point existence. The conclusion are summarized in 
Sect.IV. 



II. MODEL OF QUARK-GLUON BAGS WITH 
SURFACE TENSION 



The most convenient way to study the phase structure 
of the QGBSTM is to use the isobaric partition [241 120] 
analyzing its rightmost singularities. Hence, we assume 
that after the Laplace transform the QGBSTM2 grand 
canonical partition T, /i) generates the following iso- 
baric one: 



Z{s,T,fi) 



I 



dVe~^^ Z{V,T,ij) 



■F{s,T,fi)] 



,(1) 



where the function F{s^T^ii) includes [24] the discrete 
Fh and continuous Fq volume spectra of the bags 



F{s, T, /i) = Fh{s, T, /i) + Fq{s, T, /i) = 



(2) 



+ u{T) j% e^MT. ^)-s)v- S(T, ^)v-] (3) 

Vo 

u{T) and sq{T^ii) are continuous and, at least, double 
different iable functions of their arguments (see [24l [26] 
for details). The density of bags having mass m/., eigen 
volume v/c, baryon charge hk and degeneracy Qk is given 
by (j)k{T) Qk (l){T,mk) with 



MT) ^^Jp'dp exp 



^ gk 



27r2 



2x1/2 



(4) 



The continuous part of the volume spectrum ([3| is a gen- 
eralization of exponential mass spectrum introduced by 
Hagedorn [31] and it can be steadily derived in both the 
MIT bag model [32] and finite width QGP bag model 
[26] . The term e~^'^ describes the hard-core repulsion of 
the Van der Waals type. Z1(T, //) denotes the ratio be- 
tween the T and /i dependent surface tension coefficient 



and T (the reduced surface tension coefficient hereafter) 
which has the form 



s(r,Ai) 




^Te(m)-0, 
^Ts(/x)+0. 



(5) 



At making choice in favour of such a simple surface en- 
ergy parameterization we follow the original Fisher idea 
[5 which allows one to account for the surface energy by 
considering a mean bag of volume v and surface extent 
v^. As it has been discussed in [Ml [25] the power x < 1 
inherent in bag effective surface is a constant which, in 
principle, is different from the typical FDM and SMM 
value |. 

Let us stress here that we do not require the precise 
disappearance of I](T, /i) above the critical endpoint as 
it is usual in FDM and SMM. It was shown in [24 and 
is argued here, this point is found crucial in formulat- 
ing the statistical model with deconfining cross-over (in 
contrast with previous efforts l28 ). We would like 
also to note the negative value of the reduced surface 
tension coefficient i;(T, /i) above the T^{fi)-lme in the 
(/i, T)-plane should not be surprising. It is the well- 
known fact that in the grand canonical ensemble the sur- 
face tension coefficient includes the energy and entropy 
contributions which have the opposite signs [3 [H [33] . 
Therefore, S(T, /i) < does not mean that the surface 
energy changes the sign, but it rather signals that the 
surface entropy contribution simply exceeds the surface 
energy part and results in the negative values of surface 
free energy. In other words, the number of non- spherical 
bags of fixed volumes becomes so big that the Boltzmann 
exponent which accounts for the energy "costs" of these 
bags does not provide their suppression anymore. Such 
a situation is standard for the statistical ensembles with 
the fiuctuating extensive characteristics (the surface of 
fixed volume bag fiuctuates around the mean value) [34] . 

By construction the isobaric partition ([T]) develops two 
types of singularities: the simple pole s* = sh{T^ii) 
which is defined by the equation 



F(8*,T,/i), 



(6) 



and in addition there appears an essential singular- 
ity 5* = sq (T, ji) which is defined by the point s = 
sq (T, /i) — where the continuous part of spectrum 
FQ{s^T^jii) ([3| becomes divergent. This singularity is 
also defined by Eq.Q. Usually the statistical models 
similar to QGBSTm\24][291|30] have the following struc- 
ture of singularities. The pressure of low energy density 
phase (confined) Ph{T^ fi) is described by the simple pole 
s — sniF^n) — ^^[^'^^ which is the rightmost singu- 
larity of the isobaric partition ([T]), whereas the pressure 
of high energy density phase (deconfined) pq(T, /i) de- 
fines the system's pressure, if the essential singularity 
s = sq{T^/i) = ^^^^'^^ of this partition becomes the 
rightmost one (see Figjl]). Such an interplay of right- 
most isobaric partition singularity and the pressure of 
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the grand canonical ensemble is the typical feature of the 
Laplace transform technique [29l [30] . 

The deconfinement PT occurs at the equilibrium line 
Tc(/i) where both singularities match each other 



T = Te(M) . (7) 



In this equation one can easily recognize the Gibbs crite- 
rion for phase equilibrium. Such a behavior of the right- 
most singularities is shown in Fig.l. 

It was demonstrated in [24] the deconfinement PT 
takes place if the phase equilibrium temperature ^ is 
lower than the temperature of the null surface tension 
line ([5| for the same value of baryonic chemical potential, 
i.e. Tc(/i) < T5](/i), whereas at low values of /i the PT is 
degenerated into a cross-over because the line T = T^{/J^) 
leaves the QGP phase to appear in the hadronic phase. 
The intersection point {fiend'', Tc{jii end)) of these two lines 
Tc(/i) = Tj:{jii) is the tricritical endpoint [24l since for 
M ^ Mend and T > Tc{ fiend) at the null surface tension 
line T = Tj:{fi) there exists the surface induced PT [24 . 

The important element of our deliberation here is a 
way found out to get rid of the surface induced PT and 
to 'hide' it inside the deconfining one. In order to demon- 
strate the result we assume the surface tension coefficient 
changes its sign exactly at the deconfinement PT line, i.e. 
for max{/i(Tc)} > ja > fiend and T < T^fiend) one has 
Tc(/i) = T^ifi) while keeping the cross-over transition for 
fi < fiend similar to [24]. The possibility to match these 
two PT lines was clear long ago, but a nontriviality is seen 
in the fact that an existence of both the critical endpoint 
at {fJi end] Tc{ fiend)) and the 1*^ order deconfinement PT 
at Tc{fi) = Tj^ifi) is generated by an entire change of the 
rightmost singularity pattern. 



III. CONDITIONS FOR THE CRITICAL 
ENDPOINT EXISTENCE 

Under adopted assumption the rightmost singularity 
in the QGBSTM2 is always the simple pole since in the 
right hand side vicinity of s ^ sq (T, /i) + the value 
of Fq{s^T^ fi) ^ oo for S = S+ < 0. Then the motion 
of singularities corresponds to Fig. 2 in this situation. 
The question, however, appears whether such a behav- 
ior corresponds to PT. To clarify the point it is conve- 
nient to introduce the variable = A{Ty: ± 0, /i) 



the right most singularity s = s*{T^ — 0, /i) below and 
5+ = 5*(T5] + 0,/i) above the PT line T^fi) = T^{fi) 
for the same magnitudes of fi. Due to the relation be- 
tween the system pressure p(T, fi) and the rightmost sin- 
gularity 5*(T, fi) = ^^^T^, the difference of T derivatives, 

~T ^ ^ revealed on both sides of the PT line is 
defined by the difference of the corresponding entropy 
densities. Therefore, according to the standard classifi- 
cation of the PT order an appearance of nonzero values 



of 



a(A+-A-) 

d T 



7^ signals about the 1^^ order PT. 




FIG. 1: [Color online] Singularities of the isobaric partition 
([T]) and the corresponding graphical solution of Eq. (|6|) which 
describes a PT in the models similar to QGBSTM. The so- 
lution of Eq. ^ is shown by a filled hexagon. F{s,^) is 
shown by a solid curve for a few values of the parameter 
sets ^. F{s,^) diverges for s < sq{^) (shown by dashed 
lines), but is finite at s = sq(^) (shown by black circle). 
At low values of the parameters ^ = ^a, which can be ei- 
ther ^ = {T,fi — const} or ^ = {T = const, /j.}, the simple 
pole sh is the rightmost singularity and it corresponds to 
hadronic phase. For ^ = ^ the rightmost singular- 
ity is an essential singularity s — sq{^b), which describes 
QGP. At intermediate value ^ = both singularities coin- 
cide sh{^c) = sq{^c) and this condition is a Gibbs criterion 
([7|. At transition from the low energy density phase to the 
high density one the rightmost singularity changes from the 
simple pole to the essential singularity. 



Now using the auxiliary functions 

CO 

exp [—z + xz^] 



Ka {x) = J^^ 
VoA 



± - exp[-A±Fo-S±Fo"] 



{r-l)K 



r-l 



(8) 



(9) 



it is possible to rewrite the continuous part of volume 
spectrum ^ as FQ{s^,T,fi) = ^i(T)/^(A±, S±) inte- 
grating by parts the following integral 



- sq{T^ ± 0,/i) and to compare the T derivative of ^r(A ,11 ) 



dv 



exp [-A^v - 



Vo 



r-l 



^._i(A±,S±) 



T-V 



T- 1 



■(A±)- 



, (10) 



2 2 
X^ 

r — 1 — X 



(t-2)(t-1-x) 
/Cr-2x(a;) . 



(11) 
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FIG. 2: [Color online] Singularities of the isobaric partition 
([T]) and the corresponding graphical solution of Eq. ^ which 
describes a PT in the QGBSTM2. The legend corresponds 
to Fig. [l] In this case, however, the rightmost singularity for 
each phase is the simple pole, whereas at the PT the essen- 
tial singularity matches the simple pole due to the vanishing 
surface tension coefficient. 



Drawing Eqs. (10) and ( |ll| ) one can show the necessary 
condition of deconfinement PT existence at ^0 be- 
comes 5g(Ts,/i) = F^(5Q(TE,/i),Ts,/i) +^/(Ts)^r(0,0) 

and it provides +0, indeed. For r < 1 + 2x such a 



statement follows directly from the present form of (11), 
whereas for larger values of r exponent one needs to inte- 
grate /Ca (x)-functions in (11) while they converge at the 
lower integration limit for +0. 

With treating Eqs. (I8])-(11) one can easily find 



ds 



dT l+wJ,_i(A±,S±) 

l+w/._i(A±,S±)-^ 
which in the limit A^, gives 

.(0,0) 



(12) 



SA+ dA- 



Ulr- 



dT 



dT 



dT dT 



l + ^x/^_i(0,0) 



dFn 
ds 



(13) 



This is a remarkable result because it clearly shows in 
the present model the 1^^ order deconfinement PT does 
exist, if the T derivative of reduced surface tension coef- 
ficient has a discontinuity at the phase equilibrium line 
only! Thus, a discontinuity of the first derivative of a 
system pressure, which is a three-dimensional quantity, 
is generated by a discontinuity of the derivative of surface 
tension coefficient, which is a two-dimensional character- 
istics. In the other words, within the QGBSTM2 the de- 
confinement 1*^ order PT is just a surface induced one. 
The necessary condition for its existence is the finiteness 
of integrals /^_^(0,0) and /^_i(0,0) in ([l3|, i.e. r > 2. 

Moreover, to realize a PT from hadronic matter to 
QGP it is necessary to have at the PT line 



dA+ 

a T 



dA- 
d T 



^WT [^'q(^'M) ~ Ph{T^ fj)] > and, hence, at this line 



dT dT 



< 0. 



(14) 



Now it is clear that at the critical endpoint 
{lJiend]Tc{iJiend)) the cutropy density gap vanishes due to 
the disappearing difference — = 0. 

With the general parameterization of reduced surface 
tension coefficient which is consistent with ([5| 







(J 


rT5:(/i)-Tl 
[ Tj:ifi) \ 


r 
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s 






\T-Tj:M~ 





0, 



(15) 



■0, 



we are able to conclude about the powers and the 
values of coefficients > 0. It is obvious from (12) 



that > 1^ otherwise the corresponding entropy den- 
sity is divergent at the PT line. If, for instance, ("^ = 1, 
as predicted by the Hills and Dales model [33 , then 
C~ = 1, and according to (14) one has (T+ > a". 
If, however, C~ > 1, then from (14) it follows that 
a+C+(T - T^{ii))^^-^ > for T ^ T^{ii) + 0. The 
latter is consistent with the equality ("^ = 1. 

It can be shown that in accordance with ([t]) the in- 
equalities 

^ + ^ - 1] + 1^5.(0,0) ^ U5._(0,0)^ 

(16) 

are the sufficient conditions of the 1*^ order PT existence 



that provide (14) and guarantee the uniqueness of solu- 
tions A^ +0 on both sides of the PT line. 

The critical endpoint {iJiend^Tc^jiend)) exists, if in its 
vicinity the difference of coefficients cr^ vanishes as 



(J — (7 ^ 



d = T - TciHend) 



f^end) 

(17) 



with (end ^1- By construction in the fi — T plane d 



as defined by (17) vanishes at the tangent line to the 
PT curve at {fiend] Tc{fiend))' As one can easily see from 
either T or /i derivative of (12) any second derivative 
of the difference A+ — A~ = at the critical endpoint 

(Mend;^c(Mend)), if = C" = (end = 1 Ouly, which 

provides the 2^^ order PT available at this point. The 
higher order PT at the critical endpoint may exist for 

Cend — 2 



IV. CONCLUSION 

Here we presented new exactly solvable model, QGB- 
STM2, (or even the class of models) which develops the 
critical endpoint at (/iend; ^c(Mend))- This model natu- 
rally explains the transformation of the 1^^ order decon- 
fining PT into a weaker PT at the endpoint and into a 
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cross-over at low baryonic densities as driven by negative 
surface tension coefficient of the QGP bags at high en- 
ergy densities. It sheds new light on the QGP equation of 
state suggested in Ref.[24 where it has been shown that 
the deconfined QGP phase presents itself just a single 
infinite bag whereas the cross-over phase consists of the 
QGP bags of all possible volumes and only at very high 
pressure values this phase is presented by one large (in- 
finite) bag. The important consequence of such a prop- 
erty is that the deconfined QGP phase should be sepa- 
rated from the cross-over QGP by another PT which is 
induced by the change of surface tension coefficient sign 
of large bags. Furthermore, QGBSTM teaches us that 
for the Fisher exponent the 1*^ order deconfinement PT 
exists for 1 < r < 2 only, whereas at the endpoint there 
exists the 2'^^ order PT for | < r < 2 and this point is 
the tri-critical one. 

On the other hand the important message of QGB- 
STM2 is that a solvable model of the QCD critical end- 
point can be formulated for r > 2. Technically it is 
achieved by matching the deconfinement PT line with 
the line of vanishing surface tension coefficient Tj^{/J>) for 
M ^ Mend and T < Tc{iJiend)' This step leads to new 
strong assertion that the 1^^ order PT in QGBSTM2 



is not accompanied by change of the leading singular- 
ity type as was argued earlier in Refs. [29", Thus, 
the high density QGP phase is defined by not an essen- 
tial singularity of the isobaric partition ([T]) but its simple 
pole. Similar to QGBSTM the high density phase of this 
model is defined by the QGP crossover whereas the de- 
confined matter (an interior of single infinite bag) may 
exist at the mixed phase inherent in the deconfinement 
PT only. Besides we find also that the 1*^ order decon- 
fining PT, i.e. a discontinuity of the first derivative of a 
system pressure, which is a three-dimensional quantity, 
is generated by the discontinuity of surface tension co- 
efficient derivative, which is a two-dimensional quantity. 
Thus, we explicitly show that within the present model 
the deconfinement 1^^ order PT is the surface induced 
one. 

Another distinctive feature of these results is that for 
the first time we see the critical endpoint in the model 
with the constituents of nonzero proper volume exists not 
for r < 1 as in the SMM [H [15] and not for 1 < r < 2 
as the tricritical endpoints in the SMM and in the QG- 
BSTM [24 , but for r > 2, i.e. as in FDM [5 . Perhaps, 
this feature may be helpful to distinguish experimentally 
the QCD critical endpoint from the tri-critical one. 
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